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The multi-band nonthermal emission from the supernova 
remnant RX J1713. 7-3946 
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ABSTRACT 

Nonthermal X-rays and very high-energy (VHE) 7-rays have been detected from the 
supernova remnant (SNR) RX J1713. 7-3946, and especially the recent observations 
with the Suzaku satellite clearly reveal a spectral cutoff in the X-ray spectrum, which 
directly relates to the cutoff of the energy spectrum of the parent electrons. However, 
whether the origin of the VHE 7-rays from the SNR is hadronic or leptonic is still in 
debate. We studied the multi-band nonthermal emission from RX J1713. 7-3946 based 
on a semi-analytical approach to the nonlinear shock acceleration process by including 
the contribution of the accelerated electrons to the nonthermal radiation. The results 
show that the multi-band observations on RX J1713. 7-3946 can be well explained 
in the model with appropriate parameters and the TeV 7-rays have hadronic origin, 
i.e., they are produced via proton-proton (p-p) interactions as the relativistic protons 
accelerated at the shock collide with the ambient matter. 
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1 INTRODUCTION 

Supernova remnants (SNRs) are broadly believed to be ac- 
celeration sites of the Galactic cosmic rays. X-ray observa- 
. \ tions can provide an indication of electrons being acceler- 
. ated up to multi-TeV energies and very high-energy (VHE) 
$_( ' observations indicate that the particles would be acceler- 
5^ , ated up to hundreds of TeV or more in the SNRs. How- 
ever, the origin of the VHE 7-rays from SNRs is still un- 
certain because the TeV 7-rays from each source can usu- 
ally be explained either by hadronic models in which they 
are produced via p-p interactions or by leptonic ones in 
which they are fro m inverse Compton scattering of the par- 
ent e le ctrons (e.g., Berezhko fc Vol^l2006l: Aharonian et al.l 



11 jEnomoto et all I2OO2I) and H.E.S. S. (|Aharonian et all 
I2OO6I ). Recently. I Aharonian et al.l (|2007l ) reported the three- 
year H.E.S.S. observations on the remnant. The combined 
data significantly increase statistics and the VHE spectrum 
extends over three orders of magnitude up to energies ~ 100 
TeV. The X-ray observations on the SNR have been per- 
formed with Suzaku and a wide-band X-ray spectrum (0.4- 
40 keV) with high statistics has been shown by combin- 
ing the X-ray Imaging Spectromet er (XIS) and the Hard 
X-ray Detector (H XD) spectrum jTakahashi et all I2OO8I : 
iTanaka et aLllioOa ). The X-ray spectrum shows a clear cut- 
off shape, which allows the energy spectrum of the parent 
electrons to be more clearly investigated. 



20071: IZhang fc Fand l2007l. I2OO8I: (Yamazaki et all | 200~ 



Fang fc Zhand 120081 : iFang et all l2008l : iTanaka et al.l l200g ; 
Liu et al I2OO 8|). Detailed multi-band observations in X-rays, 
TeV 7-rays and especially in the MeV/GeV band with the 
Fermi Gamma-ray Space Telescope (FGST) satellite are im- 
portant to eliminate the uncertainty on the origin of TeV 
7-rays from SNRs. 

RX J1713.7-394 6 (G347.3-0.5) is a she ll-type SNR with 
faint radio emission (iLazeiidic et al.ll2004l) an d strong non- 
thermal X-ray emission ( Slane et al.l 1999^ . TeV 7-rays 
from the SNR shell w ere first detected by CANGAROO 
l|Muraishi et al.l I2OOOD and confirmed by CAN G ARC O- 
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In this paper, we investigate the multi-band nonthermal 
emission from RX J1713. 7-3946 in the frame of nonlinear dif- 
fusive shock acceleration mechanis m, wh ich has been studied 
exten si vely both numerica ll y (e.g ., lEllii 
19961: iBerezhko fc V61l3 Il997 



Ellison et al 



l2007h 



Baring fc Jones! 



Kang fc Jone3 |2006| : 



Malkov fc Drurvl I2OO1I: 



and semi-analytically (e.g. . 
Blasi. Amato fc Capriolil l2007l : 



Amato. Blasi fc Gabicil |2008| ). The origin of the TeV 



7-rays from the s ource had been numerically studied by 
iBerezhko fc Volkl (|2008l ) using the nonlinear kinetic theory 
of cosmic ray acceleration in SNRs. They shown that the 
recent high-energy observations with H. E.S.S. and Suzaku 
agree well with their previous study (|Berezhko fc Vol^ 
I2OO6I ) on the SNR and the hadronic origin of the TeV 
photons is more favored than the leptonic one. Moreover, 
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they argued that the non-thermal X-ray emission can 
correlate with the 7-ray one due to the correlation between 
the magnetic field amplification with the accelerated 
nuclear particl es and the associated st reaming instabilities 
(see details in IBerezhko fc Volkl I2OO8I ). Here we calculate 
the spectrum of the particles accelerate d at s h ocks using 
the semi-analytical method pr oposed in iBlasil (|2002h and 
iBlasi. Gabici fc Vannoiiil (|2005l ') and assume that the ac- 
celerated electrons have the same spectrum of the protons 
up to a maximum e nergy determined by synchrotron losses 
lEUison et all I2OO0I ). and then calculate the multi-band 
nonthermal photon spectrum from RX J1713. 7-3946. The 
results show that the multi-band spectrum of the SNR can 
be well reproduced with appropriate parameters in the 
model and the observed TeV 7-rays are produced mainly 
via hadronic intera ctions. 

Very recently, iMorlino. Amato fc Blasil (|2008l ') calcu- 
lated the high-energy spectrum for RX J1713. 7-3946 in the 
context of the nonlinear particle acceleration process at 
shocks, which is similar to that u sing in this paper. How- 
ever, o ur approach is different from lMorlino. Amato fc Blasil 
(|2008l ). Firstly, the shape of the multi-band spectrum is sen- 
sitive to the Mach number of the shock; therefore, we in- 
vestigate the nonthermal spectrum from the point of the 
Mach number and find the multi-band observations on the 
remnant can be well reproduced with appropriate param- 
eters in the model. Secondly, they treated the cut off in 
the spectrum of electrons as exp[— (£'/_Ee,max)^]; however, 
we find the Suzaku observation can still be well reproduced 
with the conventional cut off as exp[— (_B/i?e,max)] in our 
approach; Finally, the multi-band observations for the SNR 
consist well with the model results with the parameters (see 
Fig|3}, To = 10^ K, which is argued to be the temperature of 
the bubble around the SNR, and Mq = 8.0, corresponding 
to a shock speed o f ~ 3000 km s~^, wherea s a val ue of 10® 
K for To is used in iMorlino. Amato fc Blasil (|2008h . 

The structure of this paper is as follows. In §2, we briefiy 
review the model used here and show our calculation results, 
and give some discussion and conclusions in §3. 



2 THE MODEL AND RESULTS 

The pitch-angle averaged steady-state distribution of the 
protons accelerated at a shock in one dimension satisfies the 
diffusive transport equation (|Malkov fc Drurvl 1200 ll : IBlasil 
l2002l : I Amato. Blasi fc Gabicill2008D . 



Oo(p) = 



d_ 

dx 



dx 

1 du df{x,p) 
3 dx^ dp 



+ Q{x,p)^0, (1) 



where the coordinate x is directed along the shock normal 
from downstream to upstream, D is the diffusion coefficient 
and u is the fluid velocity in the shock frame, which equals 
U2 downstream {x < 0) and changes continuously upstream, 
from ui immediately upstream (a: = 0"*") of the subshock 
to uo at upstream infinity {x = -l-oo). With the assumption 
that the particles are injected at immediate upstream of the 
subshock, the source function can be written as Q{x,p) = 
Qo{p)5{x). For monoenergetic injection, Qo{p) is 



(2) 



where pinj is the injection momentum, ngas.i is the gas 
density at x = 0"*" and 77 is the fraction of particles 
injected in the acceleration process. With the injection 
recipe known as thermal leakage, 77 can be described as 
4(R,uh - l)g^e-^V37r^/^ (iBlasi. Gabici fc Vannonil 



I2OO5I : I Amato. Blasi fc Gabicill2008D . where R^ub = mi/w2 
is the compression factor at the subshock and ^ is a pa- 
rameter of the order of 2-4 describing the injection mo- 
mentum of the thermal particles in the downstream region 
Pini = CPth,2). We use ^ = 3.5 as in I Amato. Blasi fc Gabicil 
20081 ) . pth,2 = {2mpkBT2y^^ is the thermal peak momen- 
tum of the particles in the downstream fiuid with tem- 
perature T2, mp is the proton mass and fee is the Boltz- 
mann constant. Assuming the heating of the gas upstream 
is adiabatic, with the conservation condition of momentum 
fiuxes between the two sides of the subshock, we can de- 
rive the relation between the temperature of the gas far 
upstream To and T2, i.e., T2 = {-ygMo / Rtot)[Rsub/ Rtot - 
l/-Rtot + {l/'yaMS){Rtot/R,uh)''']To, where Mo is the fiuid 
Mach number far upstream, Rtot ~ uo/w2 is the total com- 
pression factor, 7g is the ratio of specific heats (7g = 5/3 for 
an ideal gas). 

RX J1713. 7-3946 has been thought as a remnant of type 
Il/Ib supernova and is evol ving in a bubble with a typica l 
temperature ^ lO'' K (e.g.. IMorlino. Amato fc Blasil [20081 '). 
In this paper, we use To = 10^ K, which corresponds to a 
sound speed of ~ 370 km s~^, to calculate the multi-band 
fiux for the SNR. 

With the assumption that the diffusion is p dependent 
and therefore particles with larger momenta move farther 
away from the shock than those with lower momenta, only 
particles with momentum ^ p can reach the point Xp and 
thus the pressure of the acceler ated particles at the point 
can be described as (|Blasill2002h 



fcR,p 



4tv 



Jp 



dp'p'^v{p')fo{p') 



(3) 



where v{p) is the velocity of particles with momentum p. 
Furthermore, the particle distribut ion function fo{p) at the 
shock can be implicitly written as (|Blasill2002l ) 



Mp) 



3Rt, 



RtotU{p) 



X exp 



dp' 3RtotU{p') 



p' RtotU{p')-l 



(4) 



where ngas,o is the gas density far upstream {x = -l-oo), Pmax 
is the maximum momentum of the accelerated particles. We 
use Pmax ~ 1.0 X 10'' TUpC in this paper, which is a typical 
value for an SNR with ambient magnetic strength ~ 100 ^G 
with an average shock speed o f ~ 5000 km s~^ for a young 
SNR with an age of ~ 1000 yr (|Yamazaki et al.ll2008l '). U{p) 
can be solved using an equation deduced from the conserva- 
tion of the mass and momentum fiuxes with the boundary 
condition ?7(pinj) = Rsuh/Rtot and f/(pmax) = 1, and then a 
value of -Rsub can be achieved b y an iterati ve procedure to 
satisfy the boundary conditions (|Blasill2002l ). 

Electrons accelerated at a shock are usually treated as 
test particles because they carry little momentum and have 
little infiuence on the shock structure. Therefore, their dis- 
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Figure 1. The resulting electron (panel a) and proton (panel b) speetra and the spectral energy distribution of synchrotron emission of the 
electrons (panel c) and p-p interactions (panel d) of the protons for RX J1713. 7-3946 with d = \ kpc, -Bsnr — 100 /^G, Pmax — 1.0 X 10^ 
mpC, Emax.e = 10 Te'V for Mq = 3, ngas,o = 6.05 cm-^, K^j, = 1.32 x 10"^ (dashed line); Mq = 5, ngas,o = 0.41 cm-^, K^^ = 1.00x10"-* 
(dotted line); Mq = 8, ngas,o = 0.13 cm ~3, K^^ = 3.82 x 10 "^ (soUd Une); and Mp = 12, ngas.o = 0.08 cm'^, = 2.44 X 10~^ 

(dash -dotted line) . The ATC A radio data l lLazendic et al ]|2004l) . Suzaku X-ray data l lTanaka et al ]|2008h . and H.E.S.S. l lAharonian et al.l 
l2007t) data from the years 2003, 2004 and 2005 are also shown. 



tribution can not be obtained by general considering mo- 
mentum and energy conservations. However, electrons and 
protons have the same acceleration rate if they have the 
same upstream diffusion length, so the electron and proton 
spectral shapes sho uld be similar at superthermal energies 
(|EUison et al.ll2000l ). As a result, the electrons have the same 
spectrum of the protons up to a maximum energy deter- 
mined by synchrotron losses, and in this paper we simply 



fe{x,p) = Kepf{x,p) exp(-_B(p)/£;max,e), 



(5) 



where E{p) is the kinetic energy of the electrons, iSmax.e 
is the cutoff energy due to the synchrotron losses, and the 
electron/proton ratio Kep is treated as a parameter. 

Assuming the accelerated particles distribute homoge- 
neously and most of the emission is from downstream of 
the shock, and using the distribution function at the shock 
to represent the particle distribution in the whole emitting 
zone, the volume-averaged emissivity for photons produced 
via p-p interactions can be written as 

da{E,Ep) 



Q{E) = 47rngas / dEpJp{Ep 



dE 



(6) 



where Ep is the proton kinetic energy, jigas is the ambi- 
ent gas number density, and Jp{Ep) = vp^ fo{p)dp/dEp is 
the volume-average proton density and v is the particles' 
velocity. We use the diffe rential cross-s e ction for photons 
d<j{E, Ep)/dE presented in lKamae et al.l (|2006h to calculate 



the hadronic 7-rays produced via p-p collisions. Finally, the 
photon flux observed at the earth can be obtained with 



F{E) 



VQ{E) 
47rd2 



(7) 



where d is the distance from the earth to the source and V 
is the average emitting volume of the source and can be esti- 
mated by y » (47r/3)f?RNR // ^tot , here J?snr is the radius of 
the SNR (|EUison et al.ll2000l '). Note that the emission from 
secondary pairs produced from p-p collisions is usually 
orders of magnitude smaller than that from the primary elec - 
trons and protons for a young SNR (jPang fc Zhand [2OO8I ) . 
thus we neglect the component from the secondary pairs. 

The distance of RX J1713. 7-3946 is uncer tain and has 
been revised somet imes dAharonian et al.ll200^ . The ASCA 
X-ray observation pCoyama et al.lll997l) and N ANTEN CO 
data (jPukui et al.ll2003l : iMoriguchi et al]l2005l ) indicate the 
SNR has a distance of ~ 1 kpc corresponding to an age of 
~ 1000 yr and we use this value in this paper. The magnetic 
field in the diffuse regions where the bulk of the synchrotron 
emission is produced is ^ 100 in order to explain the 
X-ray flux ratio betwe en the diffuse and compact regions 
(|Uchivama et al...2007. : ,Tanaka et al.ll2008l) . We use 100 fiG 
to calculate the synchrotron emission of the accelerated elec- 
trons and the cutoff energy of the electrons must be ~ 10 
TeV to well reproduce the X-ray spectrum observed with 
Suzaku. 
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Figure 2. Comparisons of the resulting spectral energy distribu- 
tions of p-p interactions for RX J1713. 7-3946 for Mq = 8, pmax = 
0.8 X 10^ rripC (dotted line); pmax = 1-0 X 10^ mpC (solid line); 
Pmax = 1.5 X 10^ mpC (dashed line); pmax = 2.0 X 10^ mpC (dash- 
dotte d line) with the H.E.S.S. data points l lAharonian et al] 
[2003). 



Figure 3. The resulting spectral energy distribution of syn- 
chrotron emission (dashed line) , inverse Compton scattering (dot- 
ted line) and p-p interactions (solid line) for RX J1713. 7-3946 
with Mo = 8, Pmax = 1.3 X 10^ rripC, ngas.o = 0.12 cm~'^, 
Kep = 3.92 X 10-5. Others are the same as FiglTI 



Fig[l] illustrates the resulting electron, proton and pho- 
ton spectra for the SNR RX J1713.7-3946 with d = 1 kpc, 

BSNR = 100 /iG, Pmax = 1.0 X lO'' TTlpC, iJmax.e = 10 TcV 

for Mo = 3, ngas,o = 6.05 cm"^ Kep = 1.32 x 10"^ (dashed 
line); Mq = 5, ngas,o = 0.41 cm"^, K^p = 1.00x10"'* (dotted 
line); Mo = 8, iigas.o = 0.13 cm-^ K^p = 3.82 x 10"^ (solid 
hne); and Mo = 12, ngas.o = 0.08 cm-^ K^p = 2.44 x 10"^ 
(dash-dotted line). For each Mach number, rigas.o is normal- 
ized to make the resulting flux at 1 TeV consistent with the 
H.E.S.S. observation, and Kep is set accordin g to the Suzaku 
resul ts at 2.6 keV. The ATCA radio data (Lazendic et ahl 
|2004|). Suzaku X-ray da ta (|Tanaka et allbOOsI ). and H.E.S.S. 
( Aharonian et alllioOTl ) data from the years 2003, 2004 and 
2005 are also s hown . Note that the radio data given by 
iLazendic et al.l (|2004h is just for the northwest part of the 
shell and the data points are multiplied by a factor of two 
to account for the ern ission for the whole remnant as in 
lAharonian et al.1 (|2006l ). The concave shape of the particle 
distribution is more obvious for a higher Mach number and 
the spectrum is harder at high momenta; importantly, the 
multi-band observed spectrum can be well reproduced with 
Mo = 8. 

The maximum energy of the accelerated electrons for 
the SNR RX J1713.7-3946 is limited due to the strong syn- 
chrotron losses in our model, so the shape of the synchrotron 
emission is nearly constant with different maximum energies 
of the protons. The influence of pmax on the spectrum of p-p 
collisions is shown in Fig[2] Obviously, the TeV observations 
can be well reproduced with Pmax around 1.0 x 10^ rripC and 
the spectrum with pmax = 2.0 x 10^ mpC relatively deviates 
with the bulk of the H.E.S.S. data points. Therefore the 
maximum energy of the protons accelerated in the SNR is 
about 100 TeV. In the following, we use Pmax = 1.3 x 10^ 
rripC to calculate the multi-band spectrum for the SNR RX 
J1713. 7-3946, and a similar value of 1.26 x 10^ rUpC was cho- 
sen to best fit the high-energy o bservation with H.E.S.S. by 
iMorlino. Amato fc BlasH (|2008h . 
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Figure 4. Clear comparisons of the model results with the ob- 
served data with Suzaku (upper panel) and H.E.S.S.. Others are 
the same as Fig[3] 



FigE] shows the comparison of the calculated results for 
Mo = 8 with the multi-band observations and Fig|4] gives a 
clear version of the comparisons of the model results with 
the Suzaku and H.E.S.S. observations. The multi-band non- 
thermal observations by ATCA, Suzaku and H.E.S.S. can be 
well explained by the model with these parameters. In cal- 
culating the inverse Compston emission of the accelerated 
electrons, we use the c anonical interstell a r valu es for the seed 
photon densities as in lAh aronian et al. I (|2006l ). i.e., 0.25 eV 
cm^'^ for the cosmic microwave background, 0.5 eV cm"^ for 
optical star light and 0.5 eV cm"^ for infrared background 
light. The contribution of the inverse Compton scattering 
to the final flux in the TeV band is orders of magnitude 
smaller than that of p-p interactions, and the contribution of 
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Bremsstrahlung is negligible even compared with the inverse 
Compton scattering. Obviously, the observed TeV 7-rays are 
from hadronic interactions as the relativistic protons collide 
with the ambient matter. 



3 DISCUSSION AND CONCLUSIONS 

RX J1713. 7-3946 is a well-observed shell-type SNR in the 
radio, X-ray and TeV 7-ray bands. Relatively strong mag- 
netic field with strength ^ 100 ^iG for the regions where 
the bulk of synchrotron emission is produced has been de- 
duced from the observation in X-rays. Especially, the obser- 
vation with Suzaku gives a wide-band spectrum with ener- 
gies from 0.4 to 40 keV, and the clear cutoff shape of the ob- 
served X-ray spectrum can constrain the model parameters 
better than before. Based on the semi-analytical approach 
to the particle acceleration process at a shock, we investi- 
gated the multi-band nonthermal spectrum for the SNR by 
taking the contribution of electrons into account. The ob- 
served multi-band nonthermal spectrum for the SNR can 
be well reproduced in the model with Kep ~ 3.92 x 10^*, 
A'Iq — 8.0, corresponding to a shock speed of ~ 3000 km 
s~^ for To = 10^ K, Pmax = 1.3 X 10^ rripC, and the maxi- 
mum energy of electrons due to the synchrotron loss is ~ 10 
TeV in order to make the model result consistent with the 
X-ray observation. The corresponding energy contained by 
the accelerated protons is ~ 3 x 10^" erg, and about 15% of 
the explosion energy is transferred to protons with a usual 
explo sion energy of ~ 2 x 10^^ erg (e.g., .Bcrczhk o fc Volkl 
I2OO6I '). The resuhs show that the TeV 7-rays observed by 
H.E.S.S. are produced predominately via p-p interactions 
and the observed X-ray spectrum with Suzaku can be well 
explained as the synchrotron emission from the accelerated 
electrons (see Fig|4|. 

With ngas,o ~ 0.12 cm^^ and To = 10^ K for a young 
SNR with an age of ~ 1000 yr, thermal X-rays from hot 
electrons should be detected from the remnant assuming 
that electrons and protons are in thermal equilibrium (e.g., 
iKatz fc Waxman.,20081 ) and the lack of thermal X-ray emis- 
sion from RX J1713. 7-3946 seems to conflict with the above 
scenario in which the observed TeV 7-rays predominately 
have hadronic origin. However, the conflict can be excluded 
by assuming the electrons and protons are not in equilib- 
rium and the elec tron temperature is significant ly smaller 
than the proton's (|Morlino. Amato fc BlasilbOOSl l. 

The possibility of leptonic origin of the TeV 7- 
rays has been investiga ted by many others. For ex ample, 
lAharonian et all (|200d ) and lYamazaki et all (|2008l ) found 
that the flat TeV emission detected by H.E.S.S. can not 
be reproduced with the one-zone model in which the elec- 
tron acceleration and gamma-ray emission take place in 
the same region unless a relatively large energy density of 
the ambient soft field is used in the calculation. Alterna- 
tively, another population of electrons is necessary to re- 
produce the multi-band observations with inverse Compton 
scattering as the origin of the TeV 7- rays with relatively 
low magnetic field st rength ~ 10 fiG (jTanaka et al.l [2OO8I : 
lYamazaki et al.ll2008h. wh i ch is m uch less than the value d e- 
duced bv lUchivama et al.l (120071 ) and lTanaka et al.l (|2008l '). 

Now, we investigate the possible leptonic origin of the 
TeV photons with the model in Sec. 2. FigO shows the re- 
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Figure 5. The resulting spectral energy distribution of syn- 
chrotron emission (dashed line), inverse Compton scattering (dot- 
ted line) and p-p interactions (solid line) for RX J1713. 7-3946 
for BsNR = 10 mG, -Bmax,e = 100 TeV, ngas.o = 0.03 cm-^, 
Kep = 1.23 x 10"^. Others are the same as FigIS] 

suiting spectrum for Bsnr = 10 fiG, -Emax,e — 100 TeV, 
JT-gas.o = 0.03 cm^^, and Kep = 1.23 x 10"'^ which is normal- 
ized according to the Suzaku observation at 2.6 keV. With 
such a weak magnetic field, as obtained before, the inverse 
Compton emission can not explain the flat TeV spectrum 
observed with H.E.S.S.; moreover, the synchrotron emission 
of electrons can not well reproduce the Suzaku spectrum in 
the X-ray band and deviates the ATCA radio data. 

In summary, the multi-wavelength spectrum of the SNR 
RX J1713. 7-3946 from radio to 7-ray bands can be well mod- 
eled with -BsNR ~ 100 i-iG in the semi-analytical nonlinear 
case. VHE 7-rays from the SNR are produced predominately 
via the tt" decay in p-p collisions and the X-ray spectrum ob- 
tained with Suzaku can be well explained as the synchrotron 
emission of the accelerated electrons. On the other hand, the 
radio. X-ray and VHE 7-ray observations can not be well 
explained in the case of leptonic origin of TeV 7-rays with 
lower magnetic held strength ~ 10 fj,G using the model in 
this paper. 
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